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Abstract 

This paper presents a proof that existence of a polynomial Lyapunov function is necessary and sufficient for exponential 
stability of sufficiently smooth nonlinear ordinary differential equations on bounded sets. The main result states that if there 
exists an n-times continuously differentiable Lyapunov function which proves exponential stability on a bounded subset of R" , 
then there exists a polynomial Lyapunov function which proves exponential stability on the same region. Such a continuous 
Lyapunov function will exist if, for example, the right-hand side of the differential equation is polynomial or at least n-times 
continuously differentiable. The proof is based on a generalization of the Weierstrass approximation theorem to differentiable 
functions in several variables. Specifically, we show how to use polynomials to approximate a differentiable function in the 
Sobolev norm W^'°° to any desired accuracy. We combine this approximation result with the second-order Taylor series 
expansion to find that polynomial Lyapunov functions can approximate continuous Lyapunov functions arbitrarily well on 
bounded sets. Our investigation is motivated by the use of polynomial optimization algorithms to construct polynomial 
Lyapunov functions. 
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1 Introduction 

o 



The Weierstrass approximation theorem was proven in 1885 [23]. This result demonstrated that real- valued polyno- 
mial functions can approximate real- valued continuous functions arbitrarily well with respect to the supremum norm 
on a compact interval. Various structural generalizations of the Weierstrass approximation theorem have focused 
on generalized mappings, as in Stone [20], and on alternate topologies, as in Krein [10]. Polynomial approximation 
of differentiable functions has been studied in numerical analysis of differential and partial differential equations. 
Results of relevance include the Bramble-Hilbert Lemma [2] and its generalization in Dupont and Scott [3], and 
Jackson's theorem [8], generalizations of which can be found in the work of Timan [21]. 



The approximation of Sobolev spaces by smooth functions has been studied in a number of contexts. In Everitt and 
Littlejohn [5], the density of polynomials of a single variable in certain Sobolev spaces was discussed. Other work 
considers weighted Sobolev spaces, as in Portilla, Quintana, Rodriguez, and Touris [15]. 

The density of infinitely continuously differentiable test functions in Sobolev spaces has been studied in the context 
of partial differential equations. Important results include the Meyers-Serrin Theorem [12], extensions of which can 
be found in Adams [1] or Evans [4]. 

Recently, the density of polynomials in the space of continuous functions has been used to motivate research on 
optimization over parameterized sets of positive polynomial functions. Relevant results include an extension of the 
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Wcicrstrass approximation theorem to polynomials subject to affinc constraints [14]. Examples of algorithms for 
optimization over cones of positive polynomials include SOSTOOLS [16], which optimizes over the cone of sums of 
squares of polynomials, and Gloptipoly [7], which optimizes over the dual space to obtain bounds on the primal 
objective. 

One application of polynomial optimization which has received particular attention has been construction of poly- 
nomial Lyapunov functions for ordinary nonlinear differential equations of the form 

x{t) = f{x{t)), 

where / : M" K". Several converse theorems have shown that local exponential stability of this system implies 
the existence of a Lyapunov function with certain properties of continuity. The reader is referred to Hahn [6] 
and Krasovskii [9] for an extensive treatment of converse theorems of Lyapunov. It is not generally known, however, 
under what conditions an exponentially stable system has a polynomial Lyapunov function. 

The main conclusion of this paper is summarized in Theorem 11, where we show that if / is n-times continuously dif- 
ferentiable, then exponentially stability of / on a ball is equivalent to the existence of a polynomial Lyapunov function 
which decreases exponentially along trajectories contained in that ball. The smoothness condition is automatically 
satisfied if / is polynomial. 

To establish this converse Lyapunov result, we prove two extensions to the Weierstrass approximation theorem. In 
Theorem 5, we show that for bounded regions, polynomials can be used to approximate continuously differentiable 
multivariate functions arbitrarily well in a variety of norms, including the Sobolev norm W^'°°. This means that for 
any continuously differentiable function and any 7 > 0, we can find a polynomial which approximates the function 
with error 7 and whose partial derivatives approximate those of the function with error 7. The proof is based on a 
construction using approximations to the partial derivatives. 

Our second extension combines a second order Taylor series expansion with the Weierstrass approximation theorem 
to find polynomials which approximate functions with a pointwise weight on the error given by 

w{x) = 

X 



These two extensions are combined into the main polynomial approximation result, which is Theorem 8. The appli- 
cation to Lyapunov functions is given in Proposition 9. Proposition 9 states that if there exists a sufficiently smooth 
continuous Lyapunov function which proves exponential stability on a bounded set, then there exists a polynomial 
Lyapunov function which proves exponential stability on the same set. In Section 5, we interpret our work as a 
converse Lyapunov theorem and briefly discuss implications for the use of polynomial optimization to prove stability 
of nonlinear ordinary differential equations. 

2 Notation and Background 

Let N" denote the set of length n vectors of non- negative natural numbers. Denote the unit cube in N" by := 
{a e N" : fti G {0,1}}. For x e R", ||3;||oo = maxj and ||.t||2 = \Jx?^x. Define the unit cube in R" by 
-B := {a; € M" : ||a;]|oo < !}• Let C(Bi) be the Banach space of scalar continuous functions defined on S C R" with 
norm 

ll/l|oo:=sup||/(x)||oo. 

For operators hi : X ^ X ^\etW;^hi : X ^ X denote the sequential composition of the hi. i.e. 

/li := /ll o /l2 o • • • O hn-\ o hn- 

i 

For a sufficiently regular function / : R" ^ R and a € N", we will often use the following shorthand to denote the 
partial derivative 
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where naturally, f /dx^ = f. For f2 c M", we define the following sets of differentiable functions. 

Cj(0) := |/ : D"f e C(0) for any a e N" such that ||a||i = ^ < i.j 

C^(0) := {/ : G C(17) for any a e N" such that ||a||oo = m&xaj < «.} 

C°°{B) = C^{B) = C^{B) is the logical extension to infinitely continuously differentiable functions. Note that in 
n-dimensions, C\{B) C Cl^{B) C C™(i?). We will occasionally refer to the Banach spaces W'''P{Q), which denote the 
standard Sobolev spaces of locally summable functions u : O — > M with weak derivatives D"u G Lp{Q) for |q:|i < fc 
and norm 

\a\i<k 

The following version of the Weierstrass approximation theorem in multiple variables comes from Timan [21]. 

Theorem 1. Suppose / : M" ^ R is continuous and G cM" is compact. Then there exists a sequence of polynomials 
which converges to f uniformly in G. 

3 Approximation of Differentiable Functions 

Most of the technical results of this paper concern a constructive method of approximating a differentiable function 
of several variables using approximations to the partial derivatives of that function. Specifically, the result states 
that if one can find polynomials which approximate the partial derivatives of a given function to accuracy 7/2"', 
then one can construct a polynomial which approximates the given function to accuracy 7, and all of whose partial 
derivatives approximate those of the given function to accuracy 7. 

The difficulty of this problem is a consequence of the fact that, for a given function, the partial derivatives of that 
function are not independent. For example, we have the identities 

D^'''fi^f{x, y, z) = [fix, y, 0) + £ D^'Mf{x, y, s) ds^ = D^''^'^^ [d^'-'''^ f{x, y, z)) , 

among others. Therefore, given approximations to the partial derivatives of a function, these approximations will, in 
general, not be the partial derivatives of any function. Then the problem becomes, for each partial derivative approx- 
imation, how to extract the information which is unique to that partial derivative in order to form an approximation 
to the original function. The following construction shows how this can be done. 

Definition 2. Let X be the space of 2"-tuples of continuous functions indexed using the 2" elements a € Z". Thus 
if functions fa S for all a € 2^", then these functions define an element of X, denoted {fa}aez'<- € X. Define 
the linear map K : X ^ Cl^{B) as 

K {{fa}aeZ'^) = ^ Gafa 
aeZn 

where Ga : C{B) ^ C^(^) is given by 



Gah = ( j h 

and where 



{gijh){xi,...,Xn) 



h{xi,. . . ,a;,_i,0,Xi+i, . . . ,a;„) j = 

J^' h{xi, . . . ,Xi-i,S,Xi+i, . . . ,Xn)ds j = 1. 
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In practice, the functions represent either the partial derivatives of a function or approximations to those partial 
derivatives. The following examples illustrate the construction. 

Example: If p = K{{qa,}aez'^)j then 

p{xi,X2) = / / 9(1,1) (si,S2)rfSirfS2 

JQ JO 

+ / 9(i,o)(si,0)dsi 
Jo 

/•X2 

+ / 9(0,1) (0,S2)rfS2 

Jo 

+ 9(0,0) (0,0). 

Notice that this structure automatically gives a way of approximating the partial derivatives of p. e.g. 

d 

—p{xi,X2) = J q{i,i){xi, S2) ds2 + q(ifi){xi,0). 

If n = 3, then 

PXi PX2 rX3 

p{xi,X2,X3) = / / / 9(1,1,1) (S1,S2, S3) rfSidS2rfS3 

Jo Jo Jo 

rxi /■X2 

+ / 9(i,i,o)(si,S2,0)dsids2 
Jo Jo 

rxi i-xs 

+ / 9(l,0,l)(si,0,S3)rfSirfS3 

Jo Jo 

nX2 rxs 

+ / 9(0,1,1)(0, 82,53) ds2C!S3 

Jo Jo 

rxi 

+ / 9(i,o,o)(si,0,0)rfsi 
Jo 

rx2 

+ I 9(0,1,0) (0, S2, 0) (is2 
Jo 

rx3 

+ / 9(0,0,1) (0,0, S3) rfs3 
Jo 

+ 9(0,0,0) (0,0,0). 



The following lemma shows that when the are the partial derivatives of a function, we recover the original function. 
The proof is simple and works by using a single integral identity to repeatedly expand the function. 

Lemma 3. For v e C^(B), if fa = D^'v for all a e Z", then K{{fc}aezn) = v. 



Proof. By assumption, v G C^(i3) and so D^v exist for all a G Z". The statement then follows from application 
of the identity v = gifiV + gi,i-^v recursively for i = l,...,n. That is, we make the sequential substitutions 
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V 1-^ 31,0^^ + V g2,ov + 52,1 gf^v , ...,vi-^ QnfiV + 9n,i^v. Then we have the following expansion. 

=gi,og2,ov + gi,o92,iD^°'^'-'°^v + 52,051, 

+ g2,igi,iD^'''''-'°'^v 
=51,052,053,01^ + 5i,o52,o53,i£'^°'°'^'°''"'°^w 

+ 5i,o52,i53,oI)(°'^'°--'°^i; + 5i,052,i53,il)(°'^'^'°--°^i' 



+ 52,o5i,i53,oi5(^'°' + 52,05i,i53,i£'(i'°'i'°'-' 



,0) 



V 



+ 92,igi,mfiD^''''-°^v + 52,i5i,i53,i£'^''''''°-'°^i^ 
=K{{D"v}c,ezn), 

as desired. □ 



The following lemma states that the linear map K is Lipschitz continuous in an appropriate sense. This means that 
a small error in the partial derivatives, fa, results in a small error of the construction K{{fa}aez^) and all of its 
partial derivatives. 

Lemma 4. Suppose p = {pa}aeZ" € X and q = {qa}aeZ" G X , then 



Proof. This proof works by noticing that the map from any function to any of the partial derivatives of 
K{{fa}aGZ") is defined by the composition of the operators gij, all of which have small gain. To see this, first 
note the following 



Then for all a,/3 e Z", 



f i = j, fc = 1 

i = j, fc = 0. 



ro 



n 9i,a. 



ai < Pi for some i 
f otherwise. 



Now, for any / e C{B), it follows from the mean value theorem that for any x G B, 



l(5i,l/)(a;)| = / f{xi,...,Xi-i,l',Xi+i,...,Xn)dl' 

Jo 

< sup \f{xi,...,Xi_i,S,Xi+i,...,Xn] 

se[-i,i] 

< llflU. 



||5i,l/l|oo < ll/llc 



Thus 

for any i. Also, it is clear that for any i, 

\{9i,of){x)\ = \f{xi, . . . ,Xi-l,0,Xi+l, . . . ,Xn)\ < 
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Therefore the gij have small gain, since Hffij/Hoo < ll/lloo for any Now since for any /? e Z", 



is the composition of gij, induction can be used to prove the the following for all a, /3 € 



a/3 



< 



Therefore 



\\D^Kp-D'^Kq\\, 



Q0 



< 



E 



9^ 



< 2" ma^ \\pa - Qo 



for any /3 € Z", as desired. 



□ 



The following theorem combines Lemmas 3 and 4 with the Weierstrass approximation theorem. It says that the 
polynomials are dense in with respect to the Sobolev norm for W^'°°, among others. 

Theorem 5. Suppose v G C^(B). Then for any e > 0, there exists a polynomial p, such that 

max II II oo < e. 

Proof. Since v G C}^{B), D"v G C{B) for all a G Z". By the Weierstrass approximation theorem, there exist 
polynomials qa such that 

Let q = {qa}aeZ" and p = Kq. Since the q^ are polynomial, p is polynomial. Let / = {D°'v}aez'^- By Lemma 3, 
V = Kf. Thus by Lemma 4, we have that 

max||Ll>-£>"z;|| = m&x\\D"Kq - D'^Kf\\ 
< 2" max||g„ -D"^;|| < e. 

□ 



Theorem 5 shows that for any continuously differentiable function, /, there exists an arbitrarily good polynomial 

approximation to the function, with error defined using the norm maxag^n ||_D"/||oo- The proof can be made con- 
structive by using the Bernstein polynomials to approximate the partial derivatives. If the partial derivatives are 
Lipschitz continuous, then this method also gives explicit bounds on the error. In practice, numerical experiments 
indicate that our constructions, at least in 2 dimensions, tend to have error roughly equivalent to the standard 
Bernstein poljoiomial approximations. 
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4 Polynomial Lyapunov Functions 



In this section, we demonstrate that polynomial Lyapunov functions can be used to approximate continuous Lya- 
punov functions. To be a Lyapunov function, a polynomial approximation must satisfy certain constraints. In par- 
ticular, if u is a Lyapunov function and p is a polynomial approximation to v, then p is also a Lyapunov function if 
it satisfies an error bound of the form 

v{x) — p{x) 



< e. 

For p to prove exponential stability, the derivatives of p and v must satisfy a similar bound. The justification for 
this form of the error bound is that the error should be everywhere bounded on a compact set, but in addition must 
decay to zero near the origin. 

The idea behind our proof of the existence of such a polynomial, p, is to combine a Taylor series approximation with 
the Weierstrass approximation theorem. Specifically, a second order Taylor series expansion about a point, xq, has 
the property that the error, or residue, R, satisfies 



R{x,xo) 







as a; xq. However, the error in the Taylor series is not guaranteed to converge uniformly over an arbitrary 
compact set as the order of the expansion increases. The Weierstrass approximation theorem, on the other hand, 
gives approximations which converge uniformly on a compact set, but in general no Weierstrass approximation 
will have the residual convergence property mentioned above for any point, Xq. Our approach, then, is to use a 
second order Taylor series expansion to guarantee accuracy near the origin. We then use a Weierstrass polynomial 
approximation to the error between the Taylor series and the function away from the origin to cancel out this error 
and guarantee a uniform bound. We then use an approach similar to that taken in Lemmas 4 and 5 to show that 
the map K can be used to construct polynomial approximations to differentiable functions in the norm 



max 



D°'v{x) 



We begin by combining the second order Taylor series expansion and the Weierstrass approximation. 
Lemma 6. Suppose v G Ci{B). Then for any e > 0, there exists a polynomial p such that 

p{x) — v{x) 



< e. 



Proof. Let the polynomial m be defined using the second order Taylor series expansion for v about a; = as 

mix) = .(0) + E-»a^(0) + 2 E 

i=l ^ij=l 

Then m approximates v near the origin and specifically 

for i,j = l,...,n. 
Now define 

''(^)-\v_M-r^ otherwise. 
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Then from Taylor's theorem(See, e.g. [11]), we have that 

n 1 " Q"^ 

v{x) =v{0) + J2 ■^^^(0) + ^ E ^^^j^:ririo) + R2{x) 



^ dxi 2 ^ dxidxj 



where ^ as a; ^ 0. Therefore 



h(x) = 



v{x) — m{x) R2{x) 



as a; — > and so h{x) is continuous at 0. Since v{x) — m{x) and x^x are continuous and x^x ^ on every domain 
not containing x = and every point x ^ has a neighborhood not containing a; = 0, we conclude that h{x) is 
continuous at every point x € M". 

We can now use the Weierstrass approximation theorem, which states that there exists some polynomial q such that 

||g - h\\oo < e- 

The Taylor and Weierstrass approximations are now combined as p{x) = m{x) + q{x)x'^x. Then p is polynomial and 



p{x) — v{x) 



m{x) + q{x)x^ x — v(x) 



m{x) — v{x) 



+ h{x) + {q{x) - h{x)) 



= \\q{x)-h{x)\\^<e 



□ 



The proof of the following lemma closely follows that of Lemma 4. However, the presence of the 1/x^x term poses 
significant technical challenges. In particular, small gain of the operators gij is no longer sufficient. We instead use 
an inductive reasoning, similar to small gain, which is described in the proof. 



Lemma 7. Let p = {pa}aez" € X and q = {qa}aez" € X. Then 



max 



Dl^Kp{x) - Dl^Kq{x) 
x^x 



< 2" max 
aez 



p. 



- qa{x) 



x'^x 



Proof. Recall that from the definition of gj^k, we have that 

d 



dxi 



9j,kf = < 



9j,k£-f iy^j 



which implies 



90 



n 

i=l 



9i,ai 



i=j,k = l , 
i=j,k = 



ai < Pi for some i 



f otherwise. 
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Now consider the term 



X 



We would like to obtain bounds on this function. For j = 1, and for any x G B, 



En 
k=l 



dt 



< sup 
•^e[-\xi\,\xi\] 

< sup 



1/(^1) • • • ) ^i— 1) ^5 • • • ; ^n) 

\f{...,x,_i,iy,x,+i,...)\ 



< 



Here the first inequality is due to the mean value theorem and that |a;i| < 1 and the second inequality follows since 
xf > v"^ for u G [— jxjl, Ixil]. Therefore, we have 



-{9i,\f){x) 



< 



-fix) 



Similarly, if j = 0, then for any x G B, 



-{9i,o.f){x) 



< 



< 







■ 1 ^n) 

















where the first inequality follows since xf > 0. Therefore, we have that for j G {0, 1} and i = 1, . . .n, 



:i9i,jf)ix) 



< 



Since the terms Ga are compositions of the gij, we can apply the above bounds inductively. Specifically, we see that 
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for any /? e Z", 



n / i^) 



< 



• < 



fix) 



Now that we have bounds on the G^, we can use the triangle inequality to deduce that for any £ Z"-, 



1 



K{p-q){x) 



(x) 



Pa (a;) - qa{x) 



< 2" max 



□ 



The following theorem gives the main approximation result of the paper. It combines Lemmas 6 and 7 to show that 
polynomials are dense in the space C""*"^ with respect to the weighted W^'°° norm with weight 1/x'^x, among others. 

Theorem 8. Suppose v is a function with partial derivatives 

D'^v e Cl{B) 

for all a G Z"-. Then for any e > 0, there exists a polynomial p, such that 



max 



< e. 



Proof. The proof is similar to that for Theorem 5. By assumption, € Cf{B) for all a € Z"^. By Lemma 6, there 
exist polynomial functions such that 



max 



ra{x) — D"v{x) 



e 

< — . 

- 2" 



Let r = {rdaez" and p = Kr. Then p is polynomial since the ret are polynomial. Let h — {D°'vJcteZ" ■ Then by 
Lemma 3, v = Kh. Therefore by Lemma 7, we have 
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max 



D"p{x) - D'^v{x) 



max 



D"Kr{x) - D^Khix) 



< 2" max 



ra{x) — D"v{x) 



as desired. 



□ 



We now conclude the section by using Theorem 8 to show that the existence of a sufficiently smooth Lyapunov 
function which proves exponential stability on a bounded set implies the existence of a polynomial Lyapunov function 
which proves exponential stability on the set. 

Proposition 9. Let Q C M" be bounded with radius r in norm ||-||oo imd f{x) be uniformly bounded on Bj. := {x € 
'■ \\x\\oo ^ f}- Suppose there exists a v : Br with D^v £ Cf{Br) for all a £ and such that 

Po\\x\\'^ < v{x) < 7o||a;||^ 
Vv{xff{x) < -5o\\xf, 

for some /3o > 0, 70 > and > and all x G Cl. Then for any /3 < /3o, 7 > 7o md 5 < So there exists a polynomial 
p such that 

p\\xr<p{x)<j\\xr 

Vp{xff{x) < -8\\xf 



for all x G fi. 



Proof. Let v{x) = v{rx) and 



lloo = sup |l/(x-)||oo- 
lklloo<r 



Choose < d < min{/3o — — 70, \^}- By Theorem 8, there exists a polynomial, p, such that for ||x||oo < 1, 



p{x) — v{x) 



< 



and 



< — 



for i = 1, . . . , n. Now let p{x) = p{x/r). Then for x G SI, \\x\\oo < r and so ||a;/r||oo < 1- Therefore we have the 
following for all a; G Q, 



p{x) = v{x) + p{x/r) — v{x/r) 

_ p{x/r) - v{x/r) 2 T 

— v[x) H -. , ,rp, j-T r X X 
[x/r)^ [x/r) 

> {Po - d)x^x 

> px^x. 
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Likewise, 



, , , , p{x/r) - v{x/r) 2 T 
[x/ry (x/r) 

< (70 + d)x'^x 

< 'yx'^x. 

Finally, 

^ \ {x/ry [x/r) J 

< ndb x'^x — 5qx^x 

< —5x^x. 

Thus the proposition holds for x G Q. □ 



A consequence of Proposition 9 is that when estimating exponential rates of decay, using polynomial Lyapunov 
functions does not result in a reduction of accuracy, i.e. if there exists a continuous Lyapunov function proving an 
exponential rate of decay with bound ao, then for any < a < ao, there exists a polynomial Lyapunov function 
which proves an exponential rate of decay with bound a. 



5 Lyapunov Stability 

Consider the system 

x{t) = f{x{t)) (1) 

where / : R" ^ R", /(O) = and x{0) = xq. We assume that there exists an r > such that for any ||a;o||oo < r, 
Equation (1) has a unique solution for all t>0. We define the solution map A : R" ^ ^([0, 00)) by 

{Ay){t)=x{t) 

for t>0, where x is the unique solution of Equation (1) with initial condition y. The following comes from Vidyasagar 
[22]. 

Theorem 10. Consider the system defined by Equation (1) and suppose that f G (K") for some integer k > 1. 
Suppose that there exist constants /x, i5, r > such that 

||(Axo)(i)||2 < Mlkolhe-^* 

for allt>0 and \\xo\\2 < r. Then there exists a Ci(R") function F : K" — > M and constants a, /3, 7,/U > such that 

a\\x\\l <V{x)<P\\x\\l 
^V{{Axm)<-j\\x\\l 

for all \\x\\2 < r. 

The following gives a converse Lyapunov result which may be taken as the main conclusion of the paper. 
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Theorem 11. Consider the system defined by Equation (1) where f G C""*" (M") . Suppose there exist constants 
fj,,6,r > such that 

\\Axo{t)\\2 < nWxohe-^* 

for allt>0 and \\xo\\2 < r. 

Then there exists a polynomial t; : M" ^ M and constants a, /?, 7, /it > such that 

a\\x\\l<v{x) <p\\x\\l 
Vv{xff{x)<-j\\x\\l 

for all \\x\\2 < r. 

Proof. We use Theorem 10 to prove the existence of a Lyapunov function y : R" ^ M with V G C"^^(R") satisfying 
the conditions on fl := {x ; l|.xl|2 < r}. Since Cj^(K."') C C"(]R"), Theorem 9 proves the existence of a polynomial 
function v which satisfies the theorem statement. □ 

An important corollary of Theorem 11 is that ordinary differential equations defined by polynomials have polynomial 
Lyapunov functions. Since polynomial optimization is typically applied to systems defined by polynomials, this means 
that the assumption of a polynomial Lyapunov function is not conservative. 

In polynomial optimization, it is common to use Positivstellensatz results to find locally positive polynomial Lya- 
punov functions in a manner similar to the 5-procedure. When the polynomial v can be assumed to be positive, 
i.e. v{x) > for all x, these conditions arc necessary and sufficient. Sec Stcnglc [19], Schmiidgcn [18], and Putinar 
[17] for strong theoretical contributions. Unfortunately, the polynomial Lyapunov functions are not positive since 
v{0) = 0, and so these conditions are no longer necessary and sufficient. However, Positivstellensatz results still allow 
us to search over polynomial Lyapunov functions in a manner which has proven very effective in practice. 

Definition 12. A polynomial, p, is sum- of- squares, if there exists a K > and polynomials gi ion i = 1, . . . ,K 

such that 

K 
i=l 

Proposition 13. Consider the system defined by Equation (1) where f is polynomial. Suppose there exists a poly- 
nomial : M" — > M, a constant e > 0, and sum-of-squares polynomials si, S2, ii, ^2 : K" — > R such that 

v{x) — si{x){r — x^ x) — S2{s) — ex^x = 

and 

-Vv{xY'f{x) - ti{x){r - x^x) - t2{x) - ex^x = 
Then there exist constants iJ.,5,r > such that 

||(ATo)(t)||2 </i||.xo||2e^'* 

for all t > and xq € Y{v,r) where Y[v,r) is the largest sublevel set of v contained in the ball ||x||^ < r. 

Sec Papachristodoulou and Prajna [13] for a proof and more details on using semidefinite programming to construct 
solutions to this polynomial optimization problem. 

6 Conclusion 

The main result of this paper is a proof that exponential stability of a sufficiently smooth nonlinear ordinary 
differential equation on a bounded region implies the existence of a polynomial Lyapunov function which decreases 
exponentially on the region. A corollary of this result is that ordinary differential equations defined by polynomials 
have polynomial Lyapunov functions. An important application of polynomial programming is the search for a 
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polynomial Lyapunov function which proves local exponential stability. Our results, therefore, tend to support 
continued research into improving polynomial optimization algorithms. 

In addition, as a byproduct of our proof, we were able to give a method for constructing polynomial approximations to 
differentiable functions. The interesting feature of this construction is the guaranteed convergence of the derivatives 
of the approximation. Another consequence of the results of this paper is that the polynomials are dense in C^(B) 
with respect to the Sobolev norm W^'^{B). 



References 

[I] R. A. Adams. Sobolev Spaces. Academic Press, 1975. 

[2] J.H. Bramble and S.R. Hilbert. Bounds for a class of linear functionals with application to Hermite interpolation. 

Numer. Math., 16:362-369, 1971. 
[3] T. Dupont and R. Scott. Polynomial approximation of functions in Sobolev spaces. Mathematics of Computation, 

34(150):441 463, April 1980. 
[4] L. C. Evans. Partial Differential Equations. American Mathematical Society, 1998. 

[5] W.N. Everitt and L. L. Littlejohn. Orthogonal polynomials and approximation in Sobolev spaces. Journal of 

Computational and Applied Mathematics, 48:69-90, 1993. 
[6] W. Hahn. Stability of Motion. Springer- Verlag, 1967. 

[7] D. Hcnrion and J.-B. Lassere. GloptiPoly: Global optimization over polynomials with Matlab and SeDuMi. In 

IEEE Conference on Decision and Control, pages 747-752, 2001. 
[8] D. Jackson. Uber genauigkeit der anndherung stetiger funktionen durch ganze rationale funktionen gegebenen 

grades und trigonometrische summengegebener ordnung. PhD thesis, Gottingen, 1911. 
[9] N. N. Krasovskii. Stability of Motion. Stanford University Press, 1963. 

[10] M. G. Krein. On a problem of extrapolation of A.N. Kolmogorov. Dokl. Akad. Nauk. SSSR, 46:306-309, 1945. 

[II] J. E. Marsden and A. J. Tromba. Vector Calculus. W. H. Feeman and Company, 3rd edition, 1988. 

[l2] N. G. Meyers and J. Serrin. H=W. Proceedings of the National Academy of Sciences, 51(6):1055-1056, June 
1964. 

[13] A. Papachristodoulou and S. Prajna. On the construction of Lyapmiov functions using the sum of squares 

decomposition. In Proceedings IEEE Conference on Decision and Control, 2002. 
[14] M. Peet and P.-A. Bliman. An extension of the weierstrass approximation theorem to linear varieties: Application 

to delay systems. In Seventh IFAC Workshop on Time-Delay Systems, 2007. 
[15] A. Portilla, Y. Quintana, J.M. Rodriguez, and E. Touris. Weierstrass' theorem in weighted Sobolev spaces with 

K derivatives. Rocky Mountain Journal of Mathematics, to appear. 
[16] Stephen Prajna, A. Papachristodoulou, P. Seller, and P. A. Parrilo. New developments in sum of squares 

optimization and SOSTOOLS. In Proceedings of the American Control Conference, 2004. 
[17] M. Putinar. Positive polynomials on compact semi-algebraic sets. Indiana Univ. Math. J., 42(3):969-984, 1993. 
[l8] C. Schmiidgen. The K-moment problem for compact semi-algebraic sets. Mathematische Annalen, 289(2): 

203-206, 1991. 

[19] G. Stengle. A nullstellensatz and a positivstellensatz in semialgebraic geometry. Mathematische Annalen, 207: 
87-97, 1973. 

[20] M. H. Stone. The generalized Weierstrass approximation theorem. Mathematics Magazine, pages 167-183 and 
237-254, 1948. 

[21] A.F. Timan. Theory of approximation of functions of a real variable. Pergamon Press Ltd., 1960. 
[22] M. Vidyasagar. Nonlinear Systems Analysis. Prentice-Hall, 1993. 

[23] K. Weierstrass. Uber die analytische darstellbarkeit sogenanntcr willkiirlichcr functionen einer reellen 
veranderlichen. Sitzungsberichte der Akademie zu Berlin, pages 633-639 and 789-805, 1885. 



14 



